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The determinations of |V u t,| from the exclusive branching ratios of B — > TV and Z? — > nlv tend 
to show a tension at the level of 3(7 [1]. On the theoretical side they depend on the lattice com- 
putation of the hadronic matrix elements and the B — > % form factor f+(q 2 )- To understand 
O | the tension, improved precision and a careful analysis of the systematics involved are necessary. 

Working towards this goal, we present preliminary lattice results of the ALPHA collaboration 
for the B — > n form factor f+(q 2 ) with Nf = 2 flavours of O (a) -improved Wilson fermions. Our 
computation uses HQET in the static limit, pion masses ranging down to ~ 250 MeV, large vol- 
umes with m K L > 4, three lattice spacings, and non-perturbative renormalization. We describe 
QO the techniques adopted to reduce the statistical noise (stochastic all-to-all with full time dilution) 

C • and the contamination from excited states (smearing for the B and the pion). We estimate the size 

of the chiral and continuum extrapolations. We discuss the impact our result could have to clarify 
the above mentioned discrepancy in the determination of |V u b|. 



DESY 12-152 
SFB/CPP- 12-67 
CERN-PH-TH/2012-267 



The 30 International Symposium on Lattice Field Theory - Lattice 2012, 
June 24-29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http : //pos. sissa. it/ 



B — > it form factor with Nf — 2 0(a) improved Wilson quarks 



Fabio Bernardoni 



1. Introduction 

The precise determination of the CKM matrix elements is a key for testing the Standard Model. 
Violations of CKM unitary or discrepancies between independent determinations of the same ma- 
trix element can provide hints for New Physics. At the time when we started our work, a tension at 
the level of 3a between two exclusive determinations of |V u b| existed, as reported e.g. in the PDG 
of 2012. These determinations use the branching ratio (BR) for the processes B — > nlv and B — > TV 
from experiment combined with the form factor f+(q 2 ) and the B decay constant /b, respectively, 
from the lattice. Also an inclusive determination, based on a perturbative expansion in a s and an 
expansion in 1 /nib, is possible. The results reported by the PDG 2012, computed before ICHEP 
2012, can be summarized as follows [1]: 

| Kb | = 0.00323(31) B^nlv 

|V ub | = 0.00510(47) B^tv (1.1) 
|V" ub | = 0.00441(34) incl. 

At ICHEP 2012 the Belle collaboration reported a new result for BR(fi — > tv) [8] based on a new 
set of data, obtained with a more sophisticated tagging of the B. This result, taken alone, would 
yield a value for |V u b| that is consistent with the exclusive determination from B — > n. However, 
more data and a careful inspection of the systematics involved are needed to draw more definitive 
conclusions. 

While the experimental precision in the differential decay rate for B — > nlv has by now reached 
good precision, B — > TV events are more difficult to reconstruct and an error of the order of 20% 
on the total branching ratio has to be expected. The situation on the theoretical side is the opposite: 
the lattice computation of a form factor is more challenging than a decay constant. 
The results for the determination of /b by the ALPHA collaboration have been presented in a poster 
at this conference [5]. The aim of the work described here is the determination of f+(q 2 ) in the 
same setup, which uses the same CLS configurations. In particular in our final results we will 
use fully non-perturbative renormalization and matching. At this conference we have presented the 
progress reached so far, putting some emphasis on the techniques employed to deal with the sizable 
finite-r effects and with the signal-to-noise ratio problems. 



2. General Setup 



Neglecting the lepton masses, the SM 
prediction for the B — > nlv decay rate is: 



dr Gl\V uh \ 2 



dq 2 



F|Vub| 
247T 3 



\Pn?\f+{q 2 



(2.1) 



where = pg — and the form factor 
f+(q 2 ) is defined through the Lorentz de- 
composition of the matrix element: 



(K(p K )\V^\B(p B ))=f + (q 2 ) 




Figure 1: Quark line diagram for the three point func- 
tion. 
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To extract this matrix element from lattice simulations, we consider the ratio: 

R(tn,t B ) = *"f — , (2.3) 

I e-9-** (Pu (x„)Pu (0) ) £ (Pa, (x b )P m (0) ) 



(7T(^)|^|B(p B ))= Hm P^b)^ 2 ^ 2 (2.4) 

r,fB,*jc->°° 

where P» and P/,/ are interpolating operators for the 71 and the P meson, respectively. Our conven- 
tions for £b and t n are illustrated in the quark line diagram, in Fig. 1. 

We treat the b quark in HQET, where a different Lorentz decomposition of (ti(p 71 )\V^\B(pb)} in 
the B rest frame is more convenient: 

(n\V°\B) = v^b/[| , 
(n\V k \B) = y^p k K f ± . 

In the static limit / + is proportional to f± : 

/+ 

We have computed the ratio eq. (2.3) on 
some of the CLS ensembles, as listed in Ta- 
ble 1, with different pion masses and three 
lattice spacings[9]. These simulations use 
two 0(tf)-improved Wilson quarks in the sea 
and relatively fine lattice spacings. Boundary 
conditions are periodic and the time extent is 
T = 2L. Finite volume effects are expected 
to contribute less than 1 % to the global sys- 
tematic error, since m n L > 4 on all ensem- 
bles. 

As the b quark cannot be treated relativisti- 
cally with present lattice spacings, we em- 
ploy HQET. This is an expansion of QCD in powers of 1 /mi, and can be used to describe processes 
for which the relevant momentum scale in the rest frame of the B-meson, in our case \pn\, satis- 
fies \p K \ <^mb. The main advantage of this approach is that HQET is renormalizable, so that the 
continuum limit can be taken and errors coming from this extrapolation can be estimated. Power 
divergences have to be removed non-perturbatively, and of course also the renormalization of the 
vector current should be done with non-perturbative precision. We here use the results of [6], 
where first the renormalization group invariant (RGI) renormalization factor was computed non- 
perturbatively. Then the matching factors only contain the physical m^-dependence through the 
ratio Mb /Ajjs where Mb is the RGI b-quark mass. The necessary matching factor, Cy has been 
computed at 3-loops in continuum perturbation theory [4]. For the operator, which enters in the 




(2.5) 



id 


L/a 


a [fm] 


m K [MeV] 


m n L 


A2 


32 


0.0755 


630 


7.7 


E4 


32 


0.0658 


580 


6.2 


E5 


32 




420 


4.7 


F6 


48 




310 


5.0 


07 


64 


0.0486 


270 


4.2 



Table 1: CLS ensembles used in this exploratory 
study. 
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computation of the form factor f+{q 2 ) in the static limit, the renormalization and matching reads: 
V& D = C v (M b /A m K% GI (go)Vr + 0{\/M b ) . (2.6) 

However, the accuracy of the perturbative expansion for B -physics has been questioned in [12]. 
We use this perturbative (3-loop) matching only in the present exploratory work, but will employ 
non-perturbative matching similar to [2] in our future analysis. We will then also include 0(a)- 
improvement for the vector current operator, while in the present setup we expect discretisation 
effects of 0(aa s ). 



3. Techniques to reduce contamination(s) from excited states and finite T 



Significant contributions to the 3-point 
function in eq. (2.3) may arise from a pion 
which propagates "around the world" in 
the r -periodic time direction of the lattice. 
Apart from a free pion state, also a n-B 
state, which we denote by \jcB\, contributes 
to these finite-r effects. This can be un- 
derstood by considering the transfer matrix 
representation for eq. (2.3) at suitable time 
separations such that the contamination from 
other excited states is negligible. This is 
equivalent to the requirements 

) > 1 , A B t B > 1 , 
(3.1) 

where A n and Ab are the energy shifts be- 
tween the ground state and the first excited 
state in the %- and in the B-channel, respec- 
tively. We then obtain: 



> 1.5 
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Figure 2: t K dependence for the considered ratio 
eq. (2.3), for E5 lattice (see Table 1). In these curves 
the interpolating operator for the B meson is always 
smeared. Circles and triangles correspond respec- 
tively to a local and a smeared one for the pion. The 
red curve is obtained by fitting the finite T effects with 
eq.(4.1). 



/?(W B )e £irW2 e MBW2 



{K\V^\B) + {0\V^\{KB])e 



-E„(T-2t x )-A,, 



where: 



y/l+ e -E x (T-2t x ) 



E K -E B - 



(3.2) 



(3.3) 



[JIB] — ^\nB] 

The finite-r effect is hence proportional to the matrix element (OjV^^B]). If we assume (also 
based on numerical evidence) that (0|V|[jtB]) and (ti\V^\B) have comparable magnitude and 
that Afefii pa 0, we conclude that fmite-r effects dominate for t % > T/2. The region where one 
expects to observe an approximate plateau is limited to t n < T/2 — k/E n where k depends on 
(OjV 1 ! [icB]}/ {k\V^\B). It would be optimal to consider lattices with T > 2L or with other boundary 
conditions to compute these 3-point functions [10]. Here, we just employ gaussian smeared wave 
functions both for the B and the n, to allow the extraction of the matrix element (7r|V|B) at 
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Figure 3: Left (Right): t K (?b) dependence of the ratio in eq. (2.3) at fixed ?b (t K ) for ensemble F6 (see 
Table 1 for detailed parameters) and smeared wave functions for both the B and the n. The red curve is the 
expected finite-T behaviour fitted according to eq. (4.1). 



relatively small time separations t n and ?b- 

Using smeared pion wave functions also helps to improve the signal-to-noise ratio which is of the 
order e ^ ( ~ Elt ^ m ^ t,t for the pion propagation. The effect of the smearing is shown in Fig. 2. In the 
future we will also consider the possibility to use non-spherical smearing for non-zero momentum 
pion wave functions [7]. 



4. Results 

To extract the form factor on each ensemble we have tried two different methods. The first is 
a linear two-parameter fit at a fixed ?b 

A+B tB e- E ^ T - 2 ^ 

The value of t% needs to be chosen small enough that the signal-to-noise ratio is satisfactory but 
large enough that the excited states effects can be neglected. The good description of the data (see 
red curve in Fig. (2)) which we achieve with this fit ansatz indicates that the finite-r effects are 
properly taken into account. 

The second is a weighted average over the plateau region, which has a relatively small extent due 
to the combined presence of excited states and large finite-r effects. The quality of our plateaux 
is shown in Fig. 3, where the t n (?b) dependence at fixed t# (t n ) is plotted for one of our lattices 
(F6, see Table l). We note that the smearing in the pion sector is very efficient and a plateau sets 
in at early t n , while in the B sector excited state effects are quite large and the plateau only starts 
when the noise is already quite sizable. This is quite different from the computation of the B -meson 
energy and decay constant from two-point functions [5]. 

Both our methods to extract the form factors give consistent results. All results presented in this 
work were obtained with the two-parameter fit method. We have also checked that results obtained 
at different time slices ?b (inside the plateau region), or with a different fit range in t % are compati- 
ble. The errors were computed using the methods described in [11]. 
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Figure 4: Fit of the (^-dependence of the form factor using the parametrization proposed in [3]. The left 
panel shows a fit using only the computation from 07 and the experimental data available. The figure from 
[1] in the right panel illustrates the experimental situation. Note that points at lower q 2 are more affected by 
0(a) and 1 /mj effects. 



5. Conclusions and Outlook 



The results collected so far do not yet allow an extrapolation to the continuum limit or to the 
physical pion mass. Moreover, the vector current is not yet 0(a) -improved beyond tree-level of 
perturbation theory and we work in the static approximation of HQET. Therefore, cutoff effects 
of order aE n a s and truncation errors of order 0(E J[ /m^) are present. We are currently working 
to include in our computation all 0(a) and 1 jm^ terms and to perform the matching with QCD at 
non-perturbative accuracy. 

The encouraging outcome of our present study is that we are able to extract a signal for the form 
factor on our ensembles with a precision of around 5 — 10%. Using a parametrization of the q 1 - 
dependence of the form factor [3], we have fitted the available experimental data together with our 
determination of f+(q 2 ) at the smallest available p n (largest q 2 ) from just one of our ensembles (07 
which is closest to the physical point), This parametrization is based on the so called z-expansion 
and very general principles, like analyticity, unitarity, and asymptotic freedom. In this fit | V u b | 
appears as a free parameter and can be determined with a precision of around 15%. 
With future more precise experimental data and including the 0(\jm\f) terms in our lattice analysis, 
we hope to reduce the uncertainty on | V u b | to the 5 — 10% level. Of course, eventually we will need 
to consider Nf > 2, but we would first like to verify with Nf = 2 simulations that the extrapolation 
to the continuum limit can be carried out within our framework of non-perturbatively renormalized 
HQET. 
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